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Employing symbolic dynamics for geodesic motion on the tesselated pseudosphere, the so-called 
Hadamard-Gutzwiller model, we construct extremely long periodic orbits without compromising 
accuracy. We establish criteria for such long orbits to behave ergodically and to yield reliable 
statistics for self-crossings and avoided crossings. Self-encounters of periodic orbits are reflected in 
certain patterns within symbol sequences, and these allow for analytic treatment of the crossing 
statistics. In particular, the distributions of crossing angles and avoided-crossing widths thus come 
out as related by analytic continuation. Moreover, the action difference for Sieber-Richter pairs of 
orbits (one orbit has a self-crossing which the other narrowly avoids and otherwise the orbits look 
very nearly the same) results to all orders in the crossing angle. These findings may be helpful for 
extending the work of Sieber and Richter towards a fuller understanding of the classical basis of 
quantum spectral fluctuations. 

I. INTRODUCTION 

Q 

U 

Billiards on surfaces of negative curvature were first investigated by J. Hadamard Q. The case of constant negative 
curvature, known as the pseudosphere, has enjoyed considerable popularity since Gutzwiller's realization Q of its 
potential as a paradigm of quantum chaos. Complete hyperbolicity, the availability of symbolic dynamics, the equality 
of the Lyapunov exponents of all periodic orbits, and the validity of Selberg's trace formula are among the attractive 
features of that system. Useful introductions can be found in Refs. • For a tesselation by maximally desymmetrized 
octagons (see below) Aurich and Steiner found the spectral fluctuations of the quantum energy spectrum faithful to 
£Nl ' the Gaussian orthogonal ensemble (GOE) of random-matrix theory (7j, as illustrated in Fig. 1 for the so-called form 
factor, the Fourier transform of the energy dependent two-point correlator of the density of levels. 

One of the urgent problems in quantum chaology is to understand the rather universal validity of random-matrix 
type spectral fluctuations for chaotic dynamical systems, nowadays known as the Bohigas-Giannoni-Schmit conjecture 
||. An important first step was made by Berry || on the basis of Gutzwiller's trace formula |[o| which expresses 
the oscillatory part of the density of energy levels as a sum over periodic orbits, d osc (E) = Re^ 7 A^e 13 " 1 ^ , with S-y 
the action and A 1 a classical stability amplitude of the periodic orbit 7. Berry realized that the double sum over 
periodic orbits in the form factor, which involves the building block A-yAy e 1 ( s -i~ s -i')/ h i n each summand, must draw an 
important contribution from the diagonal terms 7 = 7' since pairs of orbits with action differences larger than Planck's 
unit h can be expected to interfere destructively and thus to cancel in the form factor. For time reversal invariant 
systems each orbit 7 and its time reverse 7™ have equal action such that Berry's "diagonal approximation" generalizes 
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to include pairs of mutually time reversed orbits and then gives the time dependent form factor as K (r) = 2t + 
where r is the time in units of the so called Heisenberg time, which is given in terms of the mean level density d(E) 
as T H (E) = 2nhd(E). 

Recently, Sieber and Richter employed the pseudosphere in their pioneering move beyond the diagonal 

approximation; they found a one-parameter family of orbit pairs within which the action difference can be steered 
to zero. One orbit within each "Sieber-Richter pair" undergoes a small-angle self-crossing which the partner orbit 
narrowly avoids. The form factor receives the contribution K^ s (t) = — 2r 2 from the family. Together with Berry's 
fQiagfV) = 2r we thus have a semiclassical understanding of at least the first two terms in the Taylor series of the 
random-matrix form factor |l3[ Kqoe{t) = 2t — rln(l + 2r). The search is now on for further families of orbit pairs 
which might yield the higher-order terms of the expansion. At the same time, one would like and will eventually have 
to go beyond the pseudosphere, in order to find the conditions for universal behavior; interesting first steps have been 
taken for quantum graphs (lj] and other billiards Jl5[ . 

In the present paper we remain with the pseudosphere for a thorough investigation of action correlations which we 
feel necessary as a basis for further progress towards an understanding of the random-matrix like spectral fluctuations 
in this prototypical dynamical system (and beyond). We shall make extensive use of symbolic dynamics in (i) 
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establishing a certain local character of the relation between an orbit and its symbol sequence (only a section of the 
symbol sequence is necessary to determine the associated section of the orbit), (ii) constructing very long periodic 
orbits (up to hundreds of thousands of symbols) with full control of accuracy, (iii) identifying the ergodic ones among 
them and establishing a simple and instructive rederivation of Huber's exponential-proliferation law, (iv) expressing 
the angle e of a self-crossing and the closest-approach distance 6 in the partner orbit of a Sieber-Richter pair in 
terms of the Mobius-transformation matrices associated with the loops joined in a crossing, (v) revealing a useful 
analytic-continuation kinship of e and S, (vi) explicitly relating the joint density P(e, l\L) for crossing angles and loop 
lengths I in orbits of total length L to the associated density P a {5, l\L) for avoided crossings, and (vii) constructing 
an expression for the action difference of a Sieber-Richter pair valid to all orders in e. 

Some of our findings are based on (overwhelming and exceptionless) numerical evidence based on large numbers of 
long periodic orbits and thus call for mathematical substantiation. 

Even though we strictly confine ourselves to the pseudosphere we expect (and have indeed begun to check) gener- 
alizability to other systems for which symbolic dynamics is available. 
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FIG. 1. Form factor K(t) for an asymmetric octagon based on energy levels determined by Aurich and Steiner M, after 
averaging over a time window At = 0.01. Smooth line for GOE. 



II. THE HADAMARD-GUTZWILLER MODEL 



The Hadamard-Gutzwiller model [^,[7| is a two-dimensional billiard on the so called Poincare disc, i.e. the unit disc 
x 2 + y 2 = \z\ 2 < 1 endowed with the metric 

, 2 A dx 2 + dy 2 dzdz* . . 

= (l-x 2 -y 2 ) 2 = (1- \z\ 2 ) 2 ■ ( ' 

The distance d(zi, z%) between two points Z\,Z2, measured along the unique geodesic connecting them, reads 

coshfe z 2 ) = 1 + (1 _^,yJ |za|a) ■ (2-2) 

The geodesies are circles intersecting the boundary at right angles. The total length of the geodesies between its two 
crossing points with the disc boundary is infinite signaling non-compactness. Free motion on that space is completely 
hyperbolic. All trajectories (geodesies) have the same Lyapunov exponent, A = 1, and that fact makes for great 
simplifications as compared to usual billiards in flat space. The Poincare disc can be regarded as the stereographic 
projection of the pseudosphere, the surface of constant negative curvature. 

The symmetry of the Poincare disc is described by the non-compact Lorentz group SU (1,1); we shall only be 
interested in its hyperbolic elements which are matrices with the structure 



M= " £ , \a\ 2 -\f3\ 2 = 1, |Rea|>l. (2.3) 



a (3 
(3* a' 

Their action on points of the complex plane is defined as the Mobius transformation 
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z' = M(z):=-^±^. (2.4) 

p*z + a* 

Inner points of the Poincare disc are transformed to inner points, and boundary points \z\ = 1 to boundary points. 
Two boundary points zs,zjj remain invariant, zs — A4(zs) and similar for zjj; the indices signaling "stable" and 
"unstable" . Using det M. = 1 we find 

1 



± v/(a + a*) 2 -4 (2.5) 



*8,v - 

where S refers to the + sign, and U to the — sign in the case Re a > 1; if Re a < —1 the opposite sign assignment 
applies. Repeated application of the transformation M. on any point z leads to the stable fixed point, Ai k (z) — > z$ 
for k — 1, 2, 3, . . ., while iteration of the inverse map leads to the unstable fixed point, Ai~ k {z) — > Zjj. 

The geodesic passing through zg and zjj (to be called the "own" geodesic of the matrix M) is invariant with respect 
to M: if z belongs to this geodesic then so does z' = M(z). The distance d(z,M.(z)) is the same for all points of the 
geodesic and equals 

I Tr M I 

L = d(z, z') = 2arccoshJ l - . (2.6) 



The inverse of a matrix M. is obtained by the replacement a — > a* , f3 — > — (3 in ( |2.3| ). These two matrices have 

their stable and unstable points interchanged. Consequently, a matrix and its inverse shift points along their common 
"own" geodesic by the same distance but in opposite directions. 

To obtain a model with compact configuration space the Poincare disc is tesselated with tiles of equal area and 
shape. Each type of tesselation is connected with a particular discrete infinite subgroup T of SU (1,1) such that 
by acting on all inner points of a tile by some fixed matrix W =/= 1 we get all inner points of some other tile; the 
boundaries are mapped to boundaries of the same or other tiles. Using all W S T we get all tiles and restore the 
whole disc from a single initial tile. A particular tile is distinguished by containing the origin z = and is called the 
fundamental domain. In fact, different tiles are identified by identifying each point z of the Poincare disc with all its 
images W(z), and so a compact configuration space is indeed arrived at. 

In the present paper we exclusively work with octagonal tiles which are all identified with the fundamental domain. 
The opposite sides of the latter are glued together such that we arrive at a Riemann surface of genus 2, i.e. a two- 
hole doughnut. Matrices W giving rise to octagonal tiles are arbitrary products of four elementary group elements 
^0: h, h & SU(1, 1) and their inverses. For simplicity, we shall mostly consider tesselation with "regular octagons " 
(Fig. 0) which have the highest possible symmetry; the pertinent elementary matrices k = 0, 1, 2, 3, are 




h= /— - ^ ' v , v Z (2-7) 



and their inverses. The inverse of Ik is easily checked to be lk+4- Since the index k describes a phase in the off-diagonal 
matrix elements we conclude 

Ik 1 = = ^(fe+4) mod 8 = h ! (2-8) 

where we have introduced k = (k + 4) mod 8. The full set of elementary matrices is thus Ik, k — 0, ... ,7. 

The sides of the boundary of the fundamental domain can now be labeled 0, 1, . . . , 7 and constructed as follows. 
Tesselation identifies opposite sides of the regular octagon as = 4, 1 = 5, 2 = 6, and 3=7. Points on side k + 4 
are thus images of their mirror symmetric counterparts on side k; the group element responsible for such imaging 
is lk+4- In particular, the points on side solve the quadratic equation U(z) — — z* and form a circle of radius 

1/2 

(2~V 2 _ 2^!) ; which is in fact a geodesic. Side k is obtained by a rotation of side by kn/8. 
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FIG. 2. Poincare disc and its tesselation with regular octagons. 

Elementary matrices other than lk+4, lead from points on side k of the fundamental domain to boundary sides in 
the next generation of tiles, which may be called "higher Brillouin zones", using an analogy with periodic lattices. 
Obviously, there are eight Brillouin zones neighboring to the fundamental domain. In all other Brillouin zones of all 
generations, opposite sides are still identified. All Brillouin zones have the same octagonal shape and the same area 
as the fundamental domain; their visual difference in size and form is caused by the metric ( |2.1| ). 

In what follows, we shall often call the elementary matrices l/. "letters", and their products lj 1 lj 2 lj 3 . . . lj n "words". 
As a convenient shorthand for a word we shall also employ just the string of indices of the letters, i.e. W = 
lj x lj 2 lj 3 . . . l jn = (Ju 32,33, ■■■3n)- 



3 y 




y 














5 




1 7 


Z 


6 





FIG. 3. Corner points identified by the group identity (2.9); starting from I and following the arrow, the identity can be 
read off. 



As already mentioned, the identification of points z = W(z) implies gluing together opposite sides of the octagon; 
the result is a Riemann surface of genus 2. On that surface, the eight corner points of the octagon in the unit disc 
coincide. The identification of the corner points is shown in Fig. a . Starting from any corner, (for instance z in 
the figure) , taking into account that elementary matrices map opposite points of the boundaries of the fundamental 
domain onto each other, we go through all corners and come back after eight steps, and read off hlohlehhhhiz) = z. 
Due to the fact that the fixed-point equation W(z) = z for W ^ has unimodular solutions \zs,u\ = 1, we conclude 
that the matrix produc t hl phlehhhh must be the identity, since the corner point z obviously fulfills z^l. In fact, 
using the explicit form (2/7) of the Ik, one checks that the foregoing product of eight matrices is the 2x2 unit matrix, 



hhhhhhhh = (5,0,3,6,1,4,7,2) 



10 
01 



(2.9) 



Any geodesic in the tesselated unit disc can be folded into the fundamental domain where it will look like a sequence 
of disjoint circular segments starting and ending on the boundaries of the octagon. Of course, when the fundamental 



4 



tile is regarded as a surface of genus 2 the circular sections in question no longer appear as disjoint: A trajectory 
leaving the fundamental domain through one of the eight sides of the octagon and reentering on the opposite side 
appears as a smooth curve on the genus 2 surface. It is only after representing the whole unit disc by a single tile with 
opposite sides glued together (or, equivalently by a surface of genus 2) that a geodesic is capable of self-crossings. 

Consider now inertial motion along the "own" geodesic of a matrix W € T. Any point z of the geodesic and its 
image W(z) ( which we know to belong to the same geodesic) are identified by the tesselation. Geodesically moving 
from z to yV( z) we are in fact traversing a periodic orbit associated with W, and the length of that periodic orbit is 
given by ( |2.6| ); the time reversed orbit is similarly associated with the matrix W _1 . Since each W E T can be written 
as a product of elementary matrices W = lj 1 lj 2 ■ ■ ■ lj n and encoded by the symbolic word h, ■ ■ ■ l n } we have in fact 
symbolic dynamics at our disposal as a tool for investigating of periodic orbits.[] 

All matrices from an equivalence class where Z is any matrix from T, have their "own" geodesies identified 

by tesselation; there is thus one and one only periodic orbit per equivalence class in T. An infinite number of words 
pertaining to the same equivalence class and thus referring to the same periodic orbit can be obtained from one 
representative word W — (ji, J2, j3, •••ir)i this is done by cyclic permutations of the letters, by similarity transforms 
(replacing W by the longer word ZWZ TR where Z is another word, and its associate Z TR = Z^ 1 stands for the 
matrix in vers e to the one represented by Z; the superscript "TR" is read as "time reversed"), and by using the group 
identity (2J)) in all its various forms. A code of the time reversed periodic orbit is produced from the original W by 
writing its symbols ji in the opposite order and replacing each ji by ji = ji + 4 mod 8. 

Like any geodesic the periodic orbit can be folded into the fundamental domain where it will consist of a finite 
number n of disjoined circular segments. After that a distinguished n-letter word can be introduced for the orbit which 
is simply the sequence of the "landing" sides of the orbit segments. (The "launching" side is always opposite to the 
"landing" side of the preceding segment and is not admitted to the symbolic code.) This encoding contains the least 
possible number of symbols among all members of the equivalence class, and is unique up to cyclic permutations. All n 
"own" geodesies of a matrix W encoded by the distinguished word and its cyclic permutations, cross the fundamental 
domain. 




FIG. 4. Two equivalent geodesies depicting the periodic orbit (3, 6); they connect the fixed points zu, zs of lzl§ and z' v , z' s 
of hh. The primitive orbit inside the octagon (bold parts of the geodesies) has two segments. 

As an example of how to construct an orbit from its symbol sequence we consider the word W = (3,6). Its 
"own" geodesic can be found calculating the matrix product I^Iq and geodesically joining the respective fixed points 
(Eq.( [2.5| )), the curve zszjj in Fig. ^. The full orbit within the fundamental domain could be obtained by folding the 
geodesic z$zu back into that domain. However, it is easier to find the "own" geodesic of the cyclicly permuted matrix 
W = Igla and joining its fixed points z' s jj. The orbit in the fundamental domain is now given by the "inner parts" 
of the two geodesies (bold intervals of z s z\j and z' s z' v in Fig. |). The "non-primed" segment of the orbit starts on 
side 2 and ends on side 3, while the "primed" one starts on 7 and ends on 6. The time reversed orbit would have the 
launching and landing sides interchanged such that its word would be (2, 7). 

Tesselation with less symmetric octagons, also corresponding to a genus 2 Riemann surface, can be implemented 
with four elementary matrices lj € SU(1, 1) and their inverses lj, chosen such that all eight matrices obey the group 



1 Note that we take "symbol" and "letter" as synonyms; neither do we distinguish "symbol sequence", "word", "equivalence 
class of words", and "orbit", unless necessary; whenever a notational distinction seems helpful we use W,A,... for words in the 
sense of symbol sequences and VV, A, . . . for the associated Mobius-transformation matrices. 
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identity ( |2.9| ). Interestingly, completely desymmetrized genus 2 surfaces do not exist: the inversion symmetry under 
z — » — z can only be destroyed for g > 3 |jl6[ , 

III. CONSTRUCTION OF LONG PERIODIC ORBITS 

We here propose a new method for constructing very long periodic orbits. Let us consider a many-letter word 
W = (ji, J2,j3, "-in) an d its matrix lj 1 lj 2 lj 3 ---lj n - In order to explicitly construct the associated periodic orbit we can 
proceed stepwise, so as to determine each of the n circular segments separately. Each step uses only part of the word 
W (the symbol assigned to a segment and its near neighbors); the necessary length of that part is determined by the 
required accuracy and not at all by the length of W. 

In the preceding section, we discussed the connection between the fixed points zs, zjj of the Mobius transformation 
associated with W = lj 1 lj 2 ij 3 -"lj„ an d the periodic orbit. We would also like to recall that the two fixed points 
determine the circular segment associated with the first letter lj x of the word; to find the next circular segment one 
has to determine the two fixed points for the equivalent word obtained by cyclically permuting lj t to the right end. 
This is how all n circular segments of the orbit in question can be found, one after the other. 

It is convenient to first determine the stable fixed point for the first circular segment. To that end, we consider 
the sequence of matrices {W\ = lj l: W2 = Zji Zj 2 > VV3 = ■ ■ ■ } which are obtained by truncating W. For each 

matrix in the sequence we solve the quadratic fixed-point equation. The sequence of unstable fixed points thus obtained 
behaves erratically. The sequence of stable fixed points, however, converges rapidly, in fact with exponentially growing 
accuracy. The limiting point itself is none other than the stable fixed point z$ of the whole word beginning with lj t . 

The convergence of the series of stable fixed points of Wk can be ascertained analytically. To that end we first observe 
that the matrix Wk has its elem ents grow exponentially with the length L of its corresponding periodic orbit. This 
is obvious from the relation (2J3) of the length of an orbit to the trace of its matrix, |Tr Wk\ = 2 cosh(Lk/2) w e ifc / 2 



and from the unimodularity of the determinant, det Wk = ota* — (3/3* — 1. While under other circumstances such 
exponential growth gives rise to inaccuracy growing out of control, we can now rejoice in the growth working in our 



favor when determining the stable fixed point of W: For large matrix elements, Eq. (2.5) simplifies to 



When proceeding to Wk+i by including the letter lj k+1 = (ll 1 \l 2 ^j and comparing the stable fixed points Zg, z s +1 of 

Wk = (p, ^* j > ^fc+i = WfcZj fc+1 , we find for their difference, using the approximation (3.1) and the unimodularity 
of det W k , 

7 k _ k+1 _ lj2 /o n\ 

Since the elements Uj are of order unity, the difference z s — z s +1 is of order e~ Lk , i.e. small once the word Wk is 
long. Exponentially fast convergence of the sequence Zg is thus obvious. 

We can now turn to the task of determining the unstable fixed point of W. As explained in the last section, the 
time reversed orbit is obtained using the reversed Mobius transformation matrix W TR = W^ 1 , which has stable 
and unstable fixed points interchanged relative to W. Consequently, we get the unstable fixed point of the circu- 
lar segment of W associated with the letter lj x as the limit of the sequence of stable fixed points of the matrices 
{^In, ^7n ^Jn-i^Jn ^In-JT -2' ■ ■ these matrices are truncations of W TR . The sequence of stable fixed points of these 
truncations converges as rapidly as the one considered before and yields the unstable fixed point zjj for the first 
circular segment of W. 

With both zu and zs determined to the desired accuracy for the segment of the orbit corresponding to lj t , one 
repeats the procedure for the subsequent circular segment, i.e. the one corresponding to lj 2 . 

As an obvious property of our method, we find that the computational effort needed to construct a periodic orbit 
with n» 1 symbols grows only linearly with n. More importantly, we only have to do local work for local properties: 
Each of the n circular segments of the orbit inside the fundamental domain is found with a number of operations 
independent of n; for each segment one has to calculate a product of a moderate number of elementary matrices (about 
18 for a 15-digit accuracy). As we proceed segmentwise, no inaccuracy can accumulate. There is thus practically 
no limit on the length of the accessible orbits. For example, we had no difficulty in calculating a periodic orbit 
corresponding to a sequence of 10 5 randomly selected symbols. 
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It may be useful to again comment on the fact that any periodic orbit is, in principle, determined by any one of its 
circular segments within the fundamental domain. One just has to take the full geodesic running between the fixed 
points zu,zs, to which the segment belongs, and fold that geodesi c into the fundamental domain. However, such 
folding is numerically unstable and cannot be implemented with finite precision for very long orbits. 

Inasmuch as we work with fixed points of Mobius transformations, our method seems to be restricted to the 
Hadamard-Gutzwiller model. However, underlying all technicalities is a locality in the relationship between orbits 
and symbolic words, and that locality does in fact make our method generalizable, as will be explained in a separate 
publication. 

IV. PRUNING 

The locality of the word-orbit relationship makes our method a promising pruning tool. Pruning generally means 
recognizing and deleting symbolic words not corresponding to physically realizable orbits p7| . In the Hadamard- 
Gutzwiller model every word W does correspond to an allowed orbit. On the other hand, there are infinitely many 
equivalent words which must be counted only once whenever sums over orbits are to be taken, as for instance in 
Selberg's trace formula. The particular word we are interested in is simply the ordered list of "landing sides" of the 
regular octagon for the succession of circular segments of the orbit (Section [nj) . The task of selecting this distinguished 
word can be regarded as a variant of pruning. 

A first step is to discard words which can be shortened, considering that the distinguished representation contains the 
least possible number of symbols; in particular, the word must not contain pairs of symbols kk with k = (fc + 4) mod 8. 
The group identity ( |2.9| ) is yet another source of "badness" : Whenever it allows to replace a word by an equivalent 
shorter one, only the latter needs to be further scrutinized. 

Real problems are due to the fact that the group identity ( |2.9| ) allows to write some four-letter parts of words 
in reversed order, e.g. (0,5,2,7) = (7,2,5,0). The average number of such reversible 4-sequences in a word with 
n random symbols is estimated in Appendix A as 0.0043n. Each of the 4-sequences must have a uniquely defined 
direction in the distinguished word; it is not known a priori, however, which direction is the "correct" one. The 
suitable word has thus to be selected among 2 0043 " candidates. For, say, n ~ 10 5 no hope can be set on any trial- 
and-error procedure like building the orbits corresponding to the equivalent words one by one and discarding those 
not lying inside the fundamental domain. 

With our algorithm we had no problem finding the correct direction of the 4-letter sequences just mentioned. 
Whenever at some stage our method produced an arc which did not cross the fundamental domain, the letter for that 
arc turned up within a reversible 4-sequence. Upon reversing that particular 4-sequence (out of the thousands present 
in the code!) the arc always returned to the fundamental domain, and construction of the periodic orbit could be 
continued. Pruning thus becomes a local problem, which can be attacked efficiently 

V. ERGODIC PERIODIC ORBITS 

It is well known that the geodesic flow in the Hadamard-Gutzwiller model is ergodic [Q such that almost all 
trajectories cover the phase space homogeneously. Introducing limited phase-space resolution we can extend the 
notion of ergodicity to periodic orbits: Almost all sufficiently long periodic orbits cover the coarse-grained phase 
space uniformly. 

The overwhelming prevalence of ergodic periodic orbits does not imply that a random sequence of uncorrelated 
symbols must yield an ergodic orbit. In fact, practically every periodic orbit so constructed is extremely non-ergodic, 
as illustrated in Fig. ||(a) by the configuration-space density for 100000 randomly chosen symbols; the grossly non- 
uniform density thus incurred is anything but ergodic. 
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(a) (b) (c) 

FIG. 5. Coverage of fundamental domain by three periodic 100000-symbol orbits: (a) Random symbol sequence yields 
extremely inhomogeneous distribution, (b) Account of two-step correlations gives almost ergodic distribution, with small 
deviations only near octagon corners, (c) Symbol sequence imported from non-periodic trajectory gives uniform density. 

In configuration space, a long periodic orbit intersects itself many times, and the distribution of self-crossing angles 
e provides another sensitive test for ergodicity. We shall consider it in some detail because of its role in the Sieber- 
Richter theory. The angle e, defined to lie in the interval < e < 7r, is complementary to the angle between the 
velocities at the crossing; see Fig. [| The number of self-crossings with crossing angles in the interval (e, e + de) in 
periodic orbits of length L yields a density P(e\L). Since no direction is distinguished the probability that an element 
dl\ of the orbit intersects with the element dli is given by the geometric projection P(e\V) oc \dl\ x d^l oc sine. 
Comparing that prediction with the plot of P(e\L), numerically obtained for our periodic orbit with 100000 randomly 
chosen symbols, we encounter blatant disagreement, despite the enormous length of the orbit. In particular, for small 
e, the number P(e\L) of self-crossings decreases like oc e 2 / 3 rather than oc e. 




90° 180° £ 



FIG. 6. Definition of crossing angle e; strongly non-ergodic density P(e\L) for orbit with 10 s random symbols. 

It is important to stress that the non-ergodic patterns in Figs. ^|(a) and |^ are not due to an unlucky choice of the 
periodic orbit by the random-number generator employed to pick symbols. As we shall see presently, two precautions 
must be taken for our way towards ergodic periodic orbits. First, we should fix the geometric length L rather than 
the number n of symbols and, second, wave good-bye to the assumption of uncorrelated symbols. 



A. Number of symbols vs. orbit length 

Within an ensemble of orbits of fixed number of symbols n the orbit length L will fluctuate. Ascribing equal 
probability to all allowed symbol sequences and invoking, for large n, the central limit theorem we have the distribution 
of L as a Gaussian, with mean and variance both proportional to n, 

g{L\n) = -^=e LL ^ L , J g(L\n)dL = 1 ; (5.1) 
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d n and A are the mean length and variance per symbol. The values of these quantities are system specific; using 
an ensemble of 100000 periodic orbits with n — 100000 randomly selected symbols each for the regular octagon we 
numerically find 

d„« 2.2568, A w 0.6283. (5.2) 

The concentration of g(L\n) in the vicinity of its maximum becomes ever more pronounced as n grows. The mean 
length per symbol for a long periodic orbit obtained by throwing honest dice for its symbols will therefore almost 
certainly be very close to d n . 

Of greater interest is a different ensemble of periodic orbits, which has the length interval (L) fixed rather than the 
number of symbols n. In particular, it is the fixed-length ensemble that one has in mind when speaking about the 
overwhelming prevalence of ergodic orbits. To find the number of orbits N(L)dL in the length interval (L, L + dL) 
we need the number of different periodic orbits with n symbols, 

u(n) = - P : s , (5.3) 
n 

where the effective number of symbols p c s — 6.98 deviates only slightly from the naive estimate 8 — 1 = 7 obtained by 
excluding kk patterns; the slight deviation is due to the group identity (^j]); the factor 1/n reflects the equivalence 
of cyclically permuted symbol sequences; a detailed discussion of p e g will be presented in Appendix ||. We proceed 
to N(L) by summing over all possibilities for the length to lie in the interval (L, L + dL) as N(L) = J2 n 9(^\ n ) l/ ( n ) 
or, approximating the sum by an integral, 

N(L) = [ dn g(L\n)v{n) = [ dn , 1 P nln P°«- (I, " a »A )2 . (5.4) 
J J V27rn 3 A 

In the interesting case of large L the integral can be evaluated using the saddle-point approximation. The saddle of 
the integrand lies at 

(5.5) 
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FIG. 7. For fixed number n of symbols, the distribution of orbit lengths is Gaussian (5.1) with maximum at < L >= nd n 
However, ergodic orbits are those at L org = nd^, in ultra-left wing of Gaussian. 



The orbit density is thus found to obey the familiar exponential proliferation law 



N(L) = ie" L . 



with the growth rate 



V : 



dn~ y/d*-2Aln P . 



'cir 



(5.6) 



(5.7) 



The latter rate must coincide with the topological entropy For billiards on surfaces of constant negative curvature 
that entropy eq uals unity and (5.6) becomes Huber's law |^|. With our approximate numerical values (5.2) inserted, 
our expression (p?7|) yields 77 s» 1.00026, in nice agreement with the value r\ — 1 required by Huber's law. 
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The orbits described by the proliferation law ( |5.6| ) are known to be ergodic in their overwhelming majority. Since 
we have just seen this law to arise mostly from contributions of orbits with the number of symbols close to n max , it 
is clear that the ergodic orbits must have the mean length per symbol given by 

L 



^/dl-2A\n PcS . (5.8) 



For the regular octagon, with its values for d n ,A given in ( p.2|) we have cIl = 1.6283. Note that this latter value 
strongly deviates from the mean length per symbol d n — 2.2568 obtained from randomly chosen sequences of n 
symbols. Therefore, orbits generated using random sequences of n symbols have an exceedingly small probability to 
be ergodic. To see t his we must realize that the orbits constituting the maximum of the integrand in the distr ibut ion 



at fixed length (5.4) (known to be ergodic) are not those forming the maximum of the Gaussian distribution (5.1) at 
fixed number of symbols, but rather belong to its ultra-left wing; the difference of the locations of the maxima of the 
two distributions in question is of order L, while the widths are only of order \[L. 

The enormous difference between the two ensembles is due to the exponential growth of the total number of orbits 
with n symbols, v{n). The orbits with the length per symbol close to di, are in extremely small proportion to all 
orbits with the same number of symbols. However, they are dominant among the orbits with a given length since 
they have many more symbols. 



The fact that the maximum of the integrand in (5.4) is formed by ergodic orbits can be regarded as due to a theorem 
of Bowen's Jig ], The number of orbits within the ergodic subset (characterized by d^) grows in the limit L — > oo as 
fast as the total number of orbits. 



B. Symbol correlations in ergodic orbits 

Correlations must be accounted for in the symbol sequence if ergodic orbits are to be constructed. We propose 
to establish such correlations with the help of geometric considerations. The number dN^ of times an ergodic orbit 
intersects a line element of infinitesimal length dl, with the crossing angle in the interval (0, <fi + d(j)), obeys 

dN^ cx sin (j> dldtfi . (5-9) 

Indeed, dN$ and the distribution of self-crossing angles P(e\L) must have the same sinusoidal angle dependence; for 
self-crossings the line element dl is a stretch of the periodic orbit itself rather than an element of some fixed line in 
the billiard. The behavior in question is also similar to the "cosine law" of optics for the angular distribution of the 
intensity of light emitted di ffus ely by some surface element. 

We now apply the rule ( |5.9| ) to an element of one side of the regular octagon, say number 0. The segment of 
an ergodic orbit launched from that element succeeds the orbit segment with the symbol 4. In search is now the 
probability for the orbit segment launched from side to land on side k, i.e. the probability P^k that a certain 
symbol k ^ follows the symbol 4 in the symbolic code. An orbit segment is uniquely determined by the position of 
its initial point on, and the crossing angle with, the "launching side" 0. We thus find the transition probability by 
first integrating ( |5.9| ) over the angle range spanned by side k with respect to the launching element, and subsequently 
integrating along side 0. Fig. || explains the pertinent geometry. 




FIG. 8. Geodesies connecting point M at side with corners divide octagon into seven sectors. Orbit starting from M with 
momentum direction within fc-th sector will arrive at side k. Integrating over all such directions and all starting points M on 
side gives transition probability P^k- 
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For an irregular octagon the procedure just explained must be applied to all launching sides fc m to get the proba- 
bilities for the symbol £; m to be followed by &g n . In the regular octagon the transition probabilities depend only on 
the absolute value of A: m — /cg n which can be 1,2 or 3; they are to be compared with the value 1/7 = 0.143 . . . which 
would correspond to the random symbol sequence without correlation. 

Alternatively, symbol correlations can be read off non-periodic trajectories with random initial conditions. Both 
approaches lead to good agreement for the single-step transition probabilities Pi^j, as documented in the following 
table, 





Pa^i 


Pi^2 


Pi^3 


Pi^i 


theory 


0.6201 


0.2089 


0.1200 


0.0510 


nonper. traj. 


0.6193 


0.2091 


0.1205 


0.0512 



(5.10) 



Symbol sequences respecting the single-step correlations just explained yield periodic orbits exhibiting ergodicity to 
an impressive degree of approximation. In particular, the average length per symbol of periodic orbits in the regular 
octagon decreases from 2.257 to 1.694, i.e. a value rather close to the ergodic limit dt — 1.628. 

To further faithfulness to ergodicity we must account for higher correlations in the symbol sequence. The next step 
involves two-step transition probabilities P,-_>j_>fc, i.e. the conditional probabilities for a certain symbol to follow any 
given two-symbol sequence. For the regular octagon two-step correlations are embodied in a set of 16 numbers which 
can be calculated upon suitably refining the geometric reasoning applied above to find the four transition probabilities 
comprising the single-step correlations. By imposing two-step correlations on the symbol sequence we brought the 
average length per symbol for periodic orbits down to 1.643; the coverage of configuration space thus obtained is 
almost uniform, as required by ergodicity; see Fig. |^(b). It is worth noting that Pi^j^k differ from the product of 
the single-step transition probabilities Pi^jPj^k by up to 30%. 

An interesting alternative way to account for all correlations of importance is to read off symbol sequences of the 
desired length from a non-periodic trajectory. The symbolic word thus obtained can be used for generating an almost 
perfectly ergodic periodic orbit which will almost everywhere coincide, to a high precision, with the corresponding 
section of the non-periodic trajectory. Only a dozen or so segments of the periodic orbit, namely those corresponding 
to the beginning and the end of the word, will deviate from the non-periodic progenitor. The ergodicity of the orbit 
thus produced will be close to perfection, as witnessed by the density distribution in Fig. (g)(c). Moreover, the two 
trajectories under discussion will have very nearly the same length per symbol; in fact, considering 10 5 orbits with 
10 5 symbols each, generated from non-periodic trajectories, w e go t the average length per sym bol as 1.614. The latter 



value differs slightly from the prediction d^ = 1.628 based on ( |5.8D and the numerical data (5.2) on the ensemble of the 
orbits with random symbol sequences; we have not bothered to look for an explanation of that (minute) discrepancy. 

To finally illustrate the suitability of the mean length per symbol as an indicator of the fidelity of periodic orbits to 
ergodicity, Fig. |^ depicts the distribution dN^/dcj) of the angle cf> between landing side and arriving periodic orbit, for 
groups of orbits of different mean lengths d per symbol. Such groups were constructed from an ensemble of periodic 
orbits pertaining to a large fixed number of random uncorrelated symbols. The distribution of lengths L within such 
an ensemble is the Gaussian depicted in Fig. [?]. The approach to ergodicity with d — ► d^ is clearly visible in Fig. || 
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FIG. 9. Number dN^/dcp of segments landing near some point of side with an gle 



for two different values of mean 



length d per symbol and for non-periodic trajectories. Fidelity to ergodic sine law (5.9) increases as d approaches dt of (5 
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VI. CROSSING ANGLES, AVOIDED-CROSSING WIDTHS, AND LOOP LENGTHS 



We can now turn to the statistics of self-crossings and associated avoided crossings which is believed to be an 
important classical ingredient in the explanation of universality of fluctuations in quantum energy spectra of chaotic 
dynamics plffMlEl . 



A. Crossings 

Each crossing divides a periodic orbit of length L into two loops whose lengths l\, I2 add up to L. In what follows, 
we shall focus on the shorter one and its length / = min(/i, I2). For the statistical considerations to follow we employ 
the number of loops p(e, l\L)de dl with lengths in the interval (I, I + dl) and crossing angles in the interval (e, e + de) 
in ergodic orbits of length L. The density of loop lengths and crossing angles, p(e,l\L), be calculated analytically for 
the Hadamard-Gutzwiller model due to the fact that each loop can be continuously deformed to a loop with crossing 
angle tt, that is, to a periodic orbit 7, as indicated in Fig. [l^. Due to the equality of the Lyapunov exponents for all 
trajectories in the Hadamard-Gutzwiller model the loop length Z 7 is uniquely determined by the crossing angle e and 
the length L 7 of the periodic orbit 7 as f[T|| 



cosh 



cosh — 
2 



cosh sin - 

2 2 



(6.1) 



for L-y large that relation simplifies to 



i 7 (e) 



21n- 



(6.2) 



We immediately infer that a loop is longer than its periodic-orbit deformation 7 by an amount independent of 7 and 
the longer the smaller the crossing angle; even for the shortest orbits wi th t he length Lq = 3.06 (the value refers to 
the regular octagon) the error of the simplified relation (fTl) relative to (6.1) is less than 2%. 




FIG. 10. Deformation of loop to periodic orbit. 



Using the basic relation (6.1), Sieber expresses p(e, l\L) as a sum over periodic orbits |ll|] 

L x—v Lj 5(1 — ly(e)) sin e 

' r„,/f-iTo 

7 



p(e,l\L) 



T 

2tvA ^ r 7V /(-l + coshi 7 )(-l + cosh£ 7 + 2 cos 2 §) 



(6.3) 



here r 7 is a repetition number equaling 1 if 7 is a primitive periodic orbit, and A — Ait is the area of the octagon. 
Due to Z 7 (e) > L 7 only orbits with L 7 < I < L/2 contribute to the sum. 

As I grows, more and more periodic orbits contribute and the sum will eventually be dominated by long orbits for 



which 2coshL 7 w expL 7 3> 1 so that the denominator may be replaced by expL 7 . Taking into account (3.2) and 



applying the sum rule of Hannay and Ozorio de Almeida |2Cfl we obtain the ergodic result, 

L 



p(e,l\L) 



irA 



I > Zo(e) , 



(6.4) 



where cosh(Zo(e/2)) = cosh(Lo/2)/ sin(e/2) gives the length of the shortest possible loop creat ed by deformation of 
the shortest periodic orbit of the system with length Lq, again according to the basic relation (|6.l[). 
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Two reduced distributions can be obtained from the general expression fl6.3p . Integrating over the angle from zero 
to some e max , we obtain the length distribution p(l\e max , L) for loops with crossing angles smaller than e max in the 
periodic orbit of length L. The resulting expression is astonishing simple, 



J5^|e m ax, L) 



p(e, l\L)de 



L 



2-kA cosh 



2 I 



y ^coth^ 



(6.5) 



<)<! 



The only e max -dependence lies in the summation condition Z 7 (e max ) < I and indicates that only those periodic orbits 7 
contribute to the sum which can be deformed to a loop of length Z 7 (e max ) smaller than the given length I. Therefore, 
we obtain a staircase function for the quantity p(l\e max , L)cosh 2 1/2. In Fig |ll| (a), (c), we display p(l\e max , for 
£max = T and e max = 7r/36, whereas (b) displays the staircase obtained from (a) after multiplication with cosh 1/2; 
these plots reveal perfect agreement between the numerical results obta ined from all loops in an ensemble of 10 6 
ergodic periodic orbits with the prediction of the periodic-orbit sum (6.5). 
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FIG. 11. Joint density of loop lengths and crossing angles; analytic result ( p.q ) for e max = n and e max = n/36 indistinguishable 
from numerical finding. Angle-dependent gap visible in (a) and (c) is origin of logarithmic correction to ergodic behavior of 
self-crossing distribution. 



Note that a gap arises in the length distribution p(l\e max , L) due to the existence of a shortest periodic orbit and 
the related shortest possible loop of angle e and length lo(e). The gap is given by < I < lo(e ma , x ). In particular, in 
Fig. |ll|(a), lo(e max — n) = Lo = 3.06 is just the length of the smallest periodic orbit of the system, whereas in Fig. |ll| 
(c), Zo(emax =_Jk) = 9.41 is the length of the shortest possible loop with angle e max = 7r/36. Generally, in accord with 
the relation (6.2) between the lengths Z 7 (e) and L 7 , the peaks related to all periodic orbits experience practically the 
same shift to the right, ~ 21n(2/e); therefore, when e max is decreased the length distribution appears to shift to the 
right rather rigidly. 

When ( |6.3[ ) is integrated over all loop lengths I, we obtain the density of crossing angles P(e\L) in ergodic periodic 
orbits of length L as 
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L/2 , „^ „ isine (L 



P(e\L) = / p(e, l\L)dl = —j- I - - l cS (e) ) . (6.6) 



The term in the right hand side proportional to L 2 could be obtained immediately using the ergodic approximation 
( |6.4[ ). The actual number of crossings is sm aller due to the already mentioned gap in the loop length distribution, 
such that the effective integration interval in ( |6.6| ) is shorter than L/2. This is taken into account by the introduction 
of the effective minimal loop length Z ff(e); that latter quantity is close to, but somewhat smaller than, the gap visible 
in Fig[ll](a), (c), i.e. the length lo(e) of the loops related to the shortest periodic orbit. The discrepancy arises from 
the shortest periodic orbits whose contribution is greatly in excess of the ergodic predi ctio n. 



A periodic-orbit sum for l e s(e) is obtained by substituting the loop length density (6.3) for the integrand in (|6.6|) 



and comparing the left and right hand sid es. Recalling that all loop lengths exceed the lengths of their associated 



periodic orbits by the same quantity (see (6.2)) and employing the sum rule, we can make the summation interval 
independent of e, 



U(e) - 21n^^ + 0(e); 
sin 5 



(e) = -- Y Ll (6.7) 

2 r iy /{-2 + 2 cosh L 7 ) (-2 + 2 cosh L 7 + 4 cos 2 e/2) 



7^2 



For long orbits (large L) the last expression is L-independent. Indeed, the increment of 0(e), when L is replaced 
by L + AL, is given by 

A0«^~ y L 7 e- L -=0, (6.8) 



■<L^< 



L + AL 



which follows from the proliferation law for periodic orbits. 

The e dependence of 0(e) comes from the shortest and thus non-generic orbits; it is so weak that for most purposes 
we can set 0(e) « 0(0) and work with 

C *(0) 

l eS (e) ~21n — — , c=2e— w3.16; (6.9) 
2 sin | 

the constant c was obtained by choosing L/2 — 12 and correspondingly allowing for periodic orbits with lengths 
L 7 < 12 in 0(0). 



B. Avoided Crossings 

Inasmuch as every periodic orbit with a small-angle self-intersections has a partner orbit which is almost identical 
save for avoiding the said crossing, one would expect the statistics of crossings and avoided crossings to be the same, at 
least in the limit of small angles. In the present subsection we shall confirm that expectation, mainly by showing that 
there is a formal analytic-continuation relationship between crossings and avoided crossings. However, it is important 
to notice that two different types of avoided crossings arise which have different properties. 

We generalize the considerations of the preceding subsection by imagining the division of a periodic orbit into two 
"loops" by an avoided crossing. As shown in Fig. [l2|, two geodesies, the intra-octagon segments of which belong to a 
periodic orbit, may either cross inside the unit disc once or not at all. As before, we let e be the angle complementary 
to the angle between the velocity vectors of the two geodesies at an intersection as defined in Fig. [12] (a) . When the 
direction of velocity in one of the two geodesies is reversed, the role of the angles (e, it — e) is interchanged. Elementary 
Euclidian geometry provides the relation between the radii ri,r2 of the two geodesic circles, the distance I between 
its midpoints and the angle e as 

sin 2 - = - ± -^- 2 , 6.10 

2 2 4rir2 

where all quantities are measured in Euclidian units. The two solutions with ± correspond to the two complementary 
angles (e, 7r — e). 

In case the two segments do not cross there exists a minimal distance 8 between them. That distance is measured 



using the metric (2.1), along the geodesic orthogonal to the said segments. The two intersection points of this geodesic 



with the segments of the orbit divide the orbit into two loops, much the same as in the case of a crossing. 
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(a) 



(b) 



FIG. 12. Geodesies in unit disc may cross once (a), or form avoided crossing. Part (b) depicts antiparallel avoided crossing; 
reversion of motion in one geodesic would interchange angles as e *-* it — e in (a) and form parallel avoided crossing in (b). 



Surprisingly, the distance <5 is given by analytic continuation (e — » i<5) of equation (6.10) 



sinh — = — 



(6.11) 



where, in contrast to 5, the radii ri,T2 and the distance I are still mean t in the Euclidian sense, as in (6.10). We do 
not want to pause and give t he st raightforward but dull proof of ( 6.11 ) her e but shall actually recover both ( 6.1C ) 
and its analytic continuation (3.11) through symbolic dynamics in Sect. VI] . 

The geodesic along which the width of an avoided crossing is measured may consist of several disjoint arcs within 



the fundamental domain; we shall come back to this complication in Sect. VII C 



In what follows, it will be important to distinguish two different types of avoided crossings: We call "antiparallel" 
the avoided crossings with antiparallel velocities at the point of minimal distance as in Fig. [l^ (b) while avoided 
crossings with parallel velocities are called "parallel". We shall be interested in the respective distributions p a (S, l\L) 
and p p (S, l\L) of the loop length I and the minimal distance 6, both meant for loops inside ergodic orbits of total 
length L; clearly, these densities correspond to the distribution of loop lengths and crossing angles considered in the 
preceding subsection. 







(a) 



(b) 



(c) 



(d) 




FIG. 13. Crossings, parallel and antiparallel avoided crossings. For avoided crossings orbit is divided into loops by geodesic 
whose length 8 is width of avoided crossing. For narrow antiparallel (c) and parallel (e) avoided crossings, loop length is almost 
same as for crossing with angle e = 8 <C 1 (b) and tt — 8 (d). 

When determining the density p(e, l\L), we exploited the idea that every loop with length I and crossing angle e in 
a given orbit of length L is related to one periodic orbit 7, of which it can be viewed as a deformation; it is that orbit 
7 which fulfills cosh(L 7 /2) = cosh(Z/2) sin(e/2). In this sense, we have the association (I, e) <-* 7. 

Similarly, each loop defined by an antiparallel avoided crossing can be viewed as a deformation of one particular 

periodic orbit 7 with length L 7 , and the same is true for parallel avoided crossings. The analytic expressions for 

the relationships 7 <-> (l a ,S) and 7' <-> (l p ,S) appear as analytic continuations of (6.1) , in analogy to the above 



continuation of the simple trigonometric identity (6.10) for the crossing angle to Q6.ll ) for the closest-approach 
distance. The precise meaning of this analytic continuation will be explained in the following section, see equation 
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( |7.2| ). Here, we just mention that antiparallel avoided crossings correspond to solutions of Eq. (6.1) with e replaced 
by id, and parallel avoided crossings to solutions with e replaced by n + i5. The lengths l a j{5) and l p -y(6) of loops 
generated by avoided crossings as defined in Fig. O(a) can be related to the lengths of periodic orbits 7 as 



= 2 arcosh 



/ cosh \ 
V sinhl J 

/cosh 



l p -yi(5) = 2 arcosh ( 



sinh I 



^ cosh I 



antiparallel , 
parallel . 



(6.12) 
(6.13) 



Of course, in general 7 and 7' are different periodic orbits. For small closest-approach distances <5, the foregoing 
expressions for the loop lengths are simplified to 



2 In- 



Lryf 



Z 7 (e) - L 7 + C(e 2 ) for e = (5 << 1 
for 5 << 1 . 



(6.14) 



In part (b) of Fig. [13], we schematically depict loops and their lengths Z 7 (e) for small crossing angles (e << 1), similarly 
in part (c) for Z a 7 (<5) with narrowly avoided antiparallel crossings (5 -C 1) , in part (d) for Z 7 (e) with large angle e~i, 
and finally in part (e) for l p j(S) with narrowly avoided parallel crossing (8 <C 1). 

In analogy to the preceding section, we proceed to the marginal distributions p p / a {l\5 max , L) defined as 



pP/ a {l\S max ,L) = 



pP/ a {S,l\L)dS. 



(6.15) 



Most interesting is the limiting case of narrowly avoided crossings. In Fig. |TJ (a) and (b) , the loop length distributions 
pQl^maxi L) for crossings and p a (l\6 nia,x , L) for antiparallel avoided crossings are displayed with e max = £ max = W18. 
Practically no difference can be observed between these two distributions, as predicted by l 1 w I" in Eq. ( 6.14 ) 
for small argument; we here encounter the expected statistical correlation between crossings and avoided crossings. 
On the other hand, for parallel avoided crossings, the distribution p p (l\8 mllx , L) with ,5 max = 71-/I8 looks dramatically 
diff erent : As displayed in (c), spikes exactly at the position of periodic orbits are observed, as predicted by l p -y{S) w L 1 
in (6.14) in the limit of narrowly avoided parallel crossings. 
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FIG. 14. Loop length distributions p(e max = 7r/16|L) for self- crossings (a) and p a (<5 max = n/W\L) for antiparallel avoided 
crossings; both practically coincide. Very different looking loop length distribution p p (<5 m ax = 7r/16|L) for parallel avoided 
crossings in (c) exhibits pe aks n ear periodic orbit lengths since narrow parallel avoided crossings create loops close to periodic 
orbits; see Fig. [13] (e) and ( 3.14 ) . 
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VII. SYMBOLIC DYNAMICS FOR CROSSINGS AND AVOIDED CROSSINGS 



A. Crossings 



We here propose to reveal how self-intersections with small angles and narrowly avoided crossings are reflected in 
symbolic words. That correspondence will lead us to some exact results for the Hadamard-Gutzwiller model, as well 
as a better understanding of the analytic-continuation relationship of crossings and avoided crossings. The ideas to 
be presented may also become important for generalizations of the Sieber-Richter theory to other models such as 
billiards, quantum graphs and maps. 

In what follows, symbol sequences with brackets {...} refer to periodic orbits, whereas sequences without brackets 
refer to parts thereof. Consider a crossing with angle e dividing a periodic orbit into two loops. The symbolic word 
of the orbit will also be divided by the crossing into two parts which refer to a loop each and will be denoted by 
A = a\a,2 ■ ■ ■ an a and B = b\b 2 . . . b nb where tij, bk are symbols from to 7; see Fig. UB. The word of the whole orbit 
will be {AB} = {aia 2 ■ ■ ■ a na b\b 2 ■ ■ ■ b nb }, with the number of symbols n = n a + n.& and the length equal to the sum 
of the loop lengths, Lab = + Zs(e). The word is assumed pruned (see Section IV). 

We recall that the loops A,B can be continuously deformed to periodic orbits {A},{_B}, whose lengths La, Lb 



related to the lengths of the loops by (6.1), i.e., as cosh w ^ — cosh 4p / sin | for W = A or B; we shall also have 
t o re late the lengths of both orbits to the traces of the corresponding Mobius-transformation matrices^] according to 

act of elementary matrices 
also holds for non-regular 



(124), i.e. cosh^ = |Tr |. Here A, B are the 2x2 matrices obtained by taking the product of elementary matrices 
(2.7), in the order given by the symbol sequ ences A, B; the foregoing expression for Lyy 
octagons, with the elementary matrices ([2.7|) replaced by the pertinent generators. 



■y = AB = aia 2 a 3 bib- 2 b s .... 




FIG. 15. Schematic decomposition of self-crossing orbit in segments, one per symbol. First letter a\ of word W and its part 
A refers to one of two segments involved in crossing; similarly, first letter 61 of B refers to other crossing segment. 



Summing up the loop lengths we get 



jab 



2 arcosh 



Tr 



AB 



2 arcosh 



Tr A/2 



sin e/2 



2 arcosh 



Tr B/2 



sine/2 



and, solving for sin 2 e/2 



e _ (Tr AB)(Tr A)(Tr B) - (Tr A) 2 - (Tr Bf 
2 ~ (Tr AB) 2 - 4 



F(A,B) 



(7.1) 



(7.2) 



except that the elementary derivation yields the first term in the numerator of the second member as 
|(Tr AB) (Tr ^4)(Tr B)\. We have, somewhat frivolously, dropped the modulus operation since numerical evidence 
suggests (Tr AB)(Tr ^4)(Tr B) > whenever the loops A, B arise from a crossing; we regret not having found a 
proof, all the more so since only the form conjectured in (7.2) allows for analytic continuation to the case of avoided 
crossings (see below) . Incidentally, the foregoing symb olic-d ynamics expression for the crossing angle is the pre vio usly 
announced variant of the elementary-geometry result (6.10). For long orbits and loops we may approximate (fT^) as 



sm — 

2 



Tr^ Trff 
TiAB 



(7.3) 



Inasmuch as a crossing is connected with the loops A, B, the right hand side of (7.2) must obey < F (A, B) < 1. 
As regards the inverse statement, one must be a little careful. There can be several divisions of the word W into parts 



In contrast to the preceding sections we here found it convenient to notationally distinguish between symbol sequences 
A, B, AB, W and the associated matrices A, £>, AB, W 
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A, B leading to e xact ly the same value of F(A, B). In this situation, the orbit has one and only one crossing with the 
angle defined by ( |7.2| ) (provided < F(A,B) < 1), and only one of the divisions mentioned actually describes the 
loops associated to the crossing. 

As an important application of the foregoing relation between crossing angle s an d symbolic words we can now find 
the condition on the word for the crossing angle to be small. As obvious from (7.3) we must have TrATrB << Tr AB. 
A crude estimate of these traces can be obtained by recalling that the length of an ergodic orbit with the (pruned!) 
word W is roughly proportional to the number nw of letters in the word, Lw ~ nyy dL, where cLl is the average 
length per symbol (close to 1.6 in the regular octagon; see section V). Combining with ( |2.6| ) we find that the trace of 
the Mobius transform W grows exponentially with nw, 



Tr W 



1 n w d L /2 



(7.4) 



If A, B are by themselves pruned words we have for the traces in (Q Tr.4 - e nAdL / 2 ,TrS 
c (n A +n B )d L /2 with the trivial conclusion sin 2 | ~ 1; the crossing angle will not be small then. 
Suppose, however, that the matrices A and B possess the structure 



3 n B d L /2 TtAB 



A = Z 1 CZ^ 1 , B = Z 2 7LZ^ X 



(7.5) 



with £,TZ generic Mobius transforms while the "insertions " Z\, Z% are products of many elementary matrices. 
Obviously the insertions disappear from the traces of the matrices A and B since Tr A = TtZ-lCZ^ 1 = Tr£ and 
s imi larly for B\ they do not cancel in the trace of the product of AB, however. Denoting Z — Z^ Z% we have from 

(PD 



sin 

and can thus approximate the crossing angle as 



Tr£ TrK 
Tr CZKZ- 1 



sin — ss — ~ e 

2 4 



-n z d L 



nz — nz 1 + nz 2 

Very long insertions Z\, Z 2 will therefore yield very small crossing angles. 



(7.6) 



(7.7) 



P(e\L) 



All symbol 
sequences 




20 40 60 80 100 120 140 160 180 Angles 

FIG. 16. Self-crossing distribution P(e\L) for ergodic orbits, and for subset of ergodic orbits with patterns ZMZ TR L excluded. 
All small-angle self-crossings disappear for latter subset. 



To understand the physics behind condition ([7^) we recall from Section II that inverting a Mobius matrix means 
time reversal for the correspon ding code; and that the time reverse of the code Z = z\ . . . z n is Z TR = z n , . .z\ with 
z = (z + 4) mod 8. Condition (7.5) for the Mobius matrices therefore means 



A = ZiLZl 



TR 



B 



Z 2 RZ 2 TR 



(7. 



The code of the total orbit can be written as {LZRZ TK }, with Z = Z]- R Z 2 ■ The shorter periodic orbits associated 
with the loops will have the codes {L} and {R}- Since cyclic permutation is allowed in the code of orbits (not in 
loops!), the insertions can be brought together to form the self-canceling sequences ZiZj R in {A} and Z 2 Zj R in {£>}. 
Using the locality of the code-to-trajectory connection (see Section III) we may say, somewhat loosely, that the loop 
A is obtained by smoothly joining the periodic orbit {L} with two pieces of trajectory, one visiting the boundaries 
of the octagon according to Z\, and the other, encoded by Z™, running through the same sequence of boundaries 
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in the opposite direction. Such is precisely the behavior of a loop with a small crossing angle: The sections in the 
vicinity of the crossing nea rly coincide up to the sense of traversal. 



The loop structure (7.8) seems in fact necessary to produce a small crossing angle. Indeed, the beginning and 
ending parts of a loop being "nearly parallel" near a small-angle crossing the corresponding symbol sequences cannot 
help bei ng, locally identical, apart from time reversal. Impressive numerical evidence for the equivalence is provided 
by Fig. gjj Depicted are the densities P(e\L) of crossing angles for a large generic set of symbol sequences and the 
subset purged of sequences of the indicated type; in the latter case no crossings with small angles are present. 

B. Avoided crossings 

The symb olic dynamics of the avoided cros sings is inv estiga ted in much the same way. One has to replace the 
loop length (6.1) by the corresponding results ( 6.1 2| ) and ( 3.13| ). That replacement suggests an interpretation as the 



analytical continuation e — ► iS and e — > tt + iS for antiparallel and parallel avoided crossings, respectively. The width 
of an avoided crossing dividing the orbit into the loops A, B becomes 

sinh 2 - = -F(A, B) « — — — antiparallel , (7.9) 

2 TiAB 

cosh 2 | = +F{A B) w + Tl ^Jf parallel , (7.10) 
2 iv AB 



where the approximate equalities refer to long orbits and loops, as in ( |7.3[ ). It is obviously necessary that F(A\ B) < 
(F(A,B) > 1) for the existence of an antiparallel (resp. par allel) avoided crossing. The width of the antiparallel 
avoided crossing will be small under the same assumptions ( [7.8| ) about the loops as in the case of small-angle crossings; 
the code of the whole orbit in both cases has the structure LZRZ TR such that the distinction can be made only after 
calculating and thus determining the sign of F(A,B). 

The Sieber-Richter theory is based on the fact that each orbit with a small crossing angle e has a partner with almost 
the same length which avoids that crossing by the width S « e (see the next subsection for precise relations between 
the loop lengths, e, and 5); the partner has practically the same loops but one of them is traversed in the opposite 
sense relative to the original orbit. We can now translate these assertions into the language of s ym bolic dynamics. As 
we have just seen, a small-angle crossing divides the orbit into loops A, B with the structure ( |7.8| ); con sequently the 
partner must consist of the loops A, B TR . Accounting for long insertions near the crossing as after ( [7.8] ) we can state 
that each Sieber-Richter pair of orbits with small e and 5 must have codes LZRZ TR and LZR TR Z^^. The crossing 
angle in one of the orbits and the closest-approach distance in its partner are defined by 

sin 2 e/2 = FiC.ZUZ- 1 ), sinh 2 6/2 = —F(C~ 1 1 Z1ZZ~ 1 ). (7.11) 

Anticipating again esi{«lwe conclude F{H,ZKZ~ X ) « -F(Cr x , ZKZ' 1 ) or, with the help of Tr£ = Tr/T 1 , 
TtCZIZZ^ 1 « — TrC^ 1 ZIZZ^ 1 . The approximate equality of the a bsolute values of the two traces is clearly necessary 
for the length of the two orbits in the pair to be close due to (|2.6|); the minus sign is less obvious. 



C. Crossing/avoided-crossing partnership 

Even though we are mainly interested in pairs of orbits with small e and i5 it is interesting to further dwell on the 
crossing/avoided-crossing partnership for large angles and widths. We shall employ symbolic dynamics to rigorously 
relate the crossing angle e to the closest- approach distance 8 for such pairs and give exact expressions for the action 
difference in a pair and the distribution P a (5\L). Such results might be of importance for future investigations of 
higher-order terms in the expansion of the form factor. 

We again imagine the symbolic code of an orbit bisected as AB; the n th e partner orbit with respect to that bisection 
will have the symbolic word AB TR \ the function F(A, B) defined in ( \7.2\ ) will become F{A, B^ 1 ) for the partner. As 
already indicated previously, numerical checks revealed the following behavior of F under this replacement: (i) If 
the bisection AB has F{A,B) > 1 the partner orbit AB TR will also have F{A,B^ 1 ) > 1; (ii) F(A,B) < entails 
< F{A,B- r ) < 1; (iii) finally, if < F{A,B) < 1 then, the partner orbit has F(A,B^ 1 ) < 0. If A or B contains 
less than 5 symbols, and F very close to 1, very rare exceptions to this rule occu r. 

Considering the previously established interpretations of the three cases (|7.9| ), ( 7.10| ), and (7J2) we can say that 



partial time inversion of an orbit has the following respective consequences: (i) A parallel avoided crossing is replaced 
by another parallel avoided crossing; (ii) an antiparallel avoided crossing becomes a crossing in the partner; (iii) a 
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crossing becomes an antiparallel avoided crossing in the partner. For very short loops with at most 5 symbols and 
angles close to n very rare exceptions occur. 

We can use the last of these findings to more precisely define an antiparallel avoided crossing |^] as the object 
emerging after the partial time reversal of some orbit AB — > AB TR (or, equivalently, AB — > A TR B ) provided the 
bisection AB has a crossing according to < F(A, B) < 1 while F{A, B~ r ) < 0. As always, the word AB is assumed 



pruned in the sense of Sect. i.e. written in the standard representation such that each of its letters indicates a side 
of the fundamental domain visited by the orbit. The code AB TB - of the partner orbit with an avoided crossing will 
then, however, not always be pruned, and the partner orbit will sometimes not have all its circular sections within 
the fundamental domain. 

An example of a bisection of a pruned word for which partial time reversal yields a non-pruned word is provided 
by A = LZ, B = RZ. Then the symbolic word of the partner AB TR = LZZ TB -R TR obviously contains mutually 
canceling parts, and LR TR may become the standard encoding of an orbit in the fundamental domain. However, 
the avoided-crossing width in this case will be determined by the function F{£Z, Z~ lr R.) in whose arguments the 
insertion Z may not be dropped. This situation occurs when the geodesic determining the width of the avoided 
crossing consists of several segments when depicted inside the octagon (see the preceding section) . There is no reason 
not to take into account avoided crossings of the type just characterized. An orbit with an avoided crossing may thus 
contain fewer symbols and be appreciably shorter than its partner with a crossing. The closer the crossing angle is to 
7r the more often one meets with that situation. 

According to our definition there is a one-to-one correspondence between crossings and avoided crossings. One 
might be tempted to infer that in each long ergodic orbit the number of crossings must be equal to the number of 
avoided crossings, but that temptation must be resisted. Given a fixed length L, there are fewer crossings than avoided 
crossings simply since by replacing a crossing with an avoided crossing one gets a shorter periodic orbit. Therefore, to 
obtain all avoided crossings with a certain width 6 in the periodic orbits of length L one has to make the crossing to 
avoided-crossing replacement in all orbits of length L + AL(S) with suitable length increments AL, and longer orbits 
are more numerous due the exponential proliferation. 

To find the relation between the density P(e\L) of crossing angles and the density P a (5\L) of the avoided-crossing 
widths we write the total number of crossings in all periodic orbits with lengths in the interval [L, L + dL] as 
N{L)P{e\L)dLde with Huber's proliferation law N(L) = c L /L. On the other hand, the number of avoided crossings 
with width 5 in these orbits will equal the number of all crossings in the "parent" orbits whose lengths lie in [L + 
AL, L + AL + dL] which is N(L + AL)P{e\L + AL)dL\^\dS. We may thus write 



P a {5\L) 



e AL 



P{e\L + AL), (7.12) 



where the approximate-equality sign again signals large L. It remains to determine the length shift AL{5) and the 
crossing angle e(5), both as functions of 8. Clearly, if we only wanted to find the ergodic part of P a (S\L) from 
P eT a(e\L ) = {L 2 /2itA) sine we could, in analogy with our procedure in Sect. drop the length shift AL everywhere 
in ( [7.12] ) except in the exponential proliferation factor; actually, we must be more ambitious and go for the next-to- 
lcading order in L as well since the latter is of relevance for the form factor. 

To deter min e the functions e AL ^ and e{5) entering the foregoing expression we invoke two consequences of the 
definition ( |2.3| ) of SU{1, 1) matrices A, B. First, we note B + B~ x = lTr B where 1 is the unit matrix. Second, on 
multiplying with A and taking the trace we get 

TtAB + TtAB^ 1 = TtA TtB . (7.13) 



Using the identity (7.13) we can rewrite the expressions (7.2) for the crossing angle and (7.9) for the avoided-crossing 
width as 



2 e (Tr.4) 2 + (Tr.4) 2 - 4 - TtAB TtAB' 1 
cos 2 - - 



2 {TtAB) 2 - 4 

, 2 S {TtA) 2 + {TtA) 2 - 4 - TtAB TtAB' 1 ,„ 

cosh - = (7.14 

2 {TtAB- 1 ) 2 -A y ' 



and conclude 



3 From here on all avoided crossing considered will be antiparallel such that the qualifier "antiparallel" will be dropped. 
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COS 

cosh 



-1\2 



2 S 



(TrAB 
(TrAB) 2 - 4 



(7.15) 



Herein employing (2.6) and denoting the lengths of the periodic orbits with crossing and avoided crossing by L t and 
Ls we proceed to 



cos 2 
cosh ■ 



sinh ■ 



sinh -f 



(7.16) 



As by now familiar, slight simplifications arise in the physically interesting case of loops long enough for TrAB to 
be much larger than both Tr„4 and TrB. The rigorous symbolic-dynamics expressions (7.14) for e and 5 then simplify 
to 



Tr AB' 1 



C0S 2^-TrAB 



cosh - w — 



Tr AB 
Tr AB^ 



(7.17) 



and yield the previously announced relation between the crossing angle in the parent orbit and the avoided crossing 
width in the partner, 



cosh - 
2 



which implies the familiar e « S for small angles. Similarly, we may combine ( 7.16| ) and (7.17) to 



TrAB- 1 _ cosh ±£ 
TxAB ~ cosh^ 



(7.18) 



(7.19) 



and thus get the length (and thus action) difference for an orbit with a self-crossing and its partner avoiding that 
crossing, 



£ S 

AL as —4 In cos - = 4 In cosh - 
2 2 



(7.20) 



the latter relation generalizes the result of Sieber and Ric hter 11 12^ , AL = e 2 /2, to all orders in e. 



An interesting geometric interpretation of our relations (7.2C jbetween e,S, and AL in the framework of hyperbolic 
geometry will be presented in Appendix |a|. It may be well to once more say that the approximate-equality sign in 
these relations points to an error which is exponentially small when the loop lengths are large. 

We can now write the factors converting the density of crossing angles into the density of avoided-crossing widths 
in J7"ll), 



e AL = 



1 



cos 



4 e ' 



de 



tanh ■ 
tanh ■ 



(7.21) 



The latter identities yield e |ff| sme = sinh^ such that the ergodic angle distribution P CIg (e\L) — (L /2ttA) sine 
yields the ergodic distribution of avoided-crossing widths P^ rg (d\L) = (L 2 /2irA) sinh<5. Indulging in higher ambitions, 
we re tain all terms of the next-to-leading (i.e. first) order w.r.t. to L contributed by (1 + ^ ± )^ 1 P(c\L + AL) in 



( f7T2l ). Taking the full angle density ([[||6j) 



P(e\L) 



nA 



21n- 



2 sini 



(7.22) 



from Sect. VI and using the above identities ( 7.21 ) to check AL — 4 ln(c/2 sin |) = — 4 ln(c/2 sinh |) we get 



P a (S\L) = -^-sinh<5 
ttA 



21n- 



2 sinh ■ 



(7.23) 



As a most welcome surprise, after combining all the ostensibly unrelated factors in (7.12), including those resulting 
from Huber's exponential proliferation law, the final result is simply the analytical continuation of the crossing 
distribution to imaginary crossing angles. 
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D. Example of a Sieber-Richter pair 



In the previous subsection we have already commented on the symbolic codes for the periodic orbits in a Sieber- 
Richter pair as {AB} and {A TR B}. To further clarify the role of symbolic dynamics for crossings, avoided crossings, 
and Sieber-Richter pairs we would like to present a concrete such pair in the regular octagon. Fig. [l7| (a) displays 
a numerically found periodic orbit with a small-angle crossing. The reader is invited to read off the orbit segments, 
starting e.g. at the side 6 in the direction shown by arrows obtaining 6^3, 7^6, 2 — > 6, 2 — > 4, 0-^5, 1 — > 2 
thus closing the loop. The starting side is always opposite to the landing side of the preceding segment, with their 
numbers differing by ±4. Hence it is sufficient to list the consecutive landing sides only: {3, 6, 6, 4, 5, 2}. The sequence 
we got is the standard symbolic code of the orbit; the underlined symbols indicate the pair of segments crossing at 
the point P. 




FIG. 17. (a): Loops A = 3, 6 (full line) and B = 6, 4, 5, 2 (dashed) make up periodic orbit AB = {3, 6, 6, 4, 5, 2}. (b): Periodic 
orbit L = {3,6} (thick line) and loop A. (c): Periodic orbit R = {4,5} (thick dashed line) and loop B. (d): Schematic picture 
of orbits L,R and longer orbit AB — {LZRZ TR } with crossing at point P obtained by deformation of L and R. 

The crossing P divides the orbit into two loops shown by the full and dashed lines in Fig.[l7| (a); the symbol 
sequences of both loops A = 3, 6 and B = 6, 4, 5, 2 start with the symbol of the segment participating in the crossing. 
The periodic orbits associated with the loops and the loops themselves are shown by thick and thin lines respectively 
in Fig. |i~7j(b) for loop A and (c) for loop B. The periodic orbit {L} associated with the loop A has the same symbol 
sequence {3, 6} as the loop itself. The code of the orbit associated with the loop B is shorter than that of the loop, 
namely, {R} = {4, 5} while B = ZRZ TR with Z = 6; of course, the insertions Z fall out of the code of the loop- 
associated orbit. Inspection of part of Fig. [l7](c) helps to appreciate insertions Z in a loop like B = ZRZ Tn as due to 
the deformation of a circular segment of the periodic orbit {R} overstepping the boundary of the octagon. Speaking 
pictorially we may say that the middle of the upper segment of R was dragged upwards such that the vertex of the 
loop crossed the boundary 2 of the octagon, reappeared at the opposite boundary 6, and finally became the crossing 
point P, hence the insertions Z = 6 and Z TR = 2. On the other hand, when deforming segment 3 of L = 3,6 to the 
loop A ( Fig. |l7](c)) the octagon boundary is never reached, and therefore no new symbols appeared in the loop A 
compared with the periodic orbit {^4}. 

To obtain the Sieber-Richter partner of the orbit in Fig. |lj (a) we have to replace one of its loops, e.g. 

A = 3, 6 by its time reverse A TR = 6, 3 = 2, 7 . The partner periodic orbit thus formed {2, 7, 6, 4, 5, 2} is shown in Fig 
HH (a) (dotted line) together with the original orbit (full line) . The partner has an antiparallel avoided crossing whose 
width is approximately equal to the crossing angle in the original orbit; the relation valid to all orders in e is given by 
( [7.18 ). Note that apart from the two segments participating in the crossing (resp. avoided crossing) the two orbits 



run very close together; thus it is these two segments which contribute mainly to the difference in the orbit lengths. 
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(a) (b) 

FIG. 18. (a) Orbit of Fig. [H] with crossing at P and partner avoiding that crossing, (b) Definition of crossing angle e and 
avoided-crossing width 8. When crossing/avoided crossing geometry (b) is folded into the fundamental domain C is identified 
with C" and D with D' , and representation (a) results. 
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APPENDIX A: ACTION DIFFERENCE FOR PAIRS FROM HYPERBOLIC TRIANGLES 



Consider a triangle XYZ formed by geodesies on the pseudosphere, with the angles X, Y, Z named after the 
respective vertices and the lengths of the respective opposite sides x, y, z. The following generalizations of elementary 
Euclidian-geometry relations hold true || 



cosh a; 



cos Y cos Z + cos X sin X sin Y 



sin Y sin Z 



sin Z 



sinh x sinh y sinh ; 



(Al) 



Let us apply these identities to the triangle FPD in Fig. 18 (b) depicting the crossing and avoided crossing in 
a Sieber-Richter pair. Again naming the angles of the triangle after the respective vertices we have F = ir/2,P = 
ir/2 — e/2. The point D is supposed to be removed to infinity (recall the assumption of long orbits) such that the lengths 
FD = L 2 and PD — Li tend to infinity while the angle D tends to zero. However, the difference L 2 — L\ remains 
finite and gives, up to an obvio us fa ctor 4, the length (or action) difference in search, AL = L e — L$ — 4(Li — L 2 ). 
From the first of the relations (Al) we get the connection between the crossing angle e and the closest-approach 
distance 5 = 2FP, 



e(S) 



1 



cosh | 



(A2) 



which we recognize as ( 7.1S ); similarly, the second group of relations in (Al) yields the length difference 



e smhLi 
cos - = . , T 
2 sinh L-2 



L1—L2 



AL sa 4 In cos • 



(A3) 



in agreement with (7.20). The association of symbolic dynamics and Mobius-transformation matrices employed in 
Sect. VII to derive the foregoing relations obviously incorporates the hyperbolic geometry of the Hadamard-Gutzwiller 
model. 



APPENDIX B: EFFECTIVE NUMBER OF SYMBOLS 



We here determine the number v(n) of non-equivalent n-letter words (alias periodic orbits with n segments) 7 with 
an alphabet containing p = 8 different letters. As mentioned in Section V, a rough estimate is v(n) w (p— l) n /n, since 
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a symbol j may not be succeeded by its "inverse" j = (j +4) mod 8, and since a factor 1/n is due to the identification 
of cyclic permutations. The influence of the group identity 

(0,5,2,7,4,1,6,3) = 1. (Bl) 

requires more thought. 

Our firs t ta sk is to find the shortest possible version of a given word. Toward that end we note that whenever the 



sequence (Bl) or any of its cyclic permutations is encountered within a word, we must delete that sequence and thus 
shorten the word by 8 letters. Similarly, if some part of the identity with m symbols and 4 < m < 8, is encountered, 
that part must be replaced by its complement with 8 — m symbols. In any such instance the word is shortened 
by 8,6,4 or 2 letters. Next, we have to consider the 4-symbol sequences for which the group identity provides an 
equivalent 4-symbol sequence like for example (0, 5, 2, 7) = (7, 2, 5, 0). If a word contains k such 4-letter patches there 
are 2 k different representations of the same orbit; only one of these representations has the property that each symbol 
corresponds to one segment of the orbit inside the fundamental domain. 

Let us start with the reversible 4-letter sequences and denote by a the probability that a symbol of a long word 
is the beginning of such a sequence. A combinatorial estimate of a can be made as follows . T here are 16 different 



4-letter patterns of the type in question, due to 8 cyclic permutations of both the identity in (Bl) and its inverse, and 
each of these patterns occurs with probability 1/(8 x 7 3 ). However, the pattern can only be included when the first 
and last letter are not the inverse of the neighboring one in the complete word, such that a "junction factor" (7/8) 2 
must be incorporated. The estimate in search thus comes out as a com b = (7/8) 2 x 16/(8 x 7 3 ) « 0.00446. Assuming 
independent occurrence of the 4-letter sequences we expect to find exactly k of them in an n-letter word with the 
binomial probability 

P n (k,a) = jz^ct{l-ar-*. (B2) 

We refrain from attempts at improving that combinatorial estimate by accounting for the "interaction" of the 4-letter 
patches caused by their finite length. Upon checking many long wor ds with randomly chosen letters we found such 



interaction effects quite unimportant: The binomial distribution (B2) is borne out very well with a num = 0.00425. 

Similarly, let (3 be the probability to find, in any letter, the beginning of any of the previously discussed m-symbol 
patterns with 4 < m < 8 which has an equivalent shorter partner with 8 — m symbols. Again assuming independence 
of several such events we get a binomial distribution P n (k, (3) as before and the combinatorial estimate for (3 is, again 
neglecting "interactions" , Aheor = (7/8) 2 [16/(7 x 8 7 ) + 16/(7 x 8 6 ) + 16/(7 x 8 5 ) + 16/(7 x 8 4 )] = 0.00074. Here again, 
the interactions were checked to be unimportant by going through a large sample of random words. The numerical 
data for P™ m (fc) so obtained once more fit the binomial distribution, with /3 n um = 0.000716. . 

We can finally put all pieces together. The number of n-letter words allowed after excluding jj sequences of symbols 
and accounting for the equivalence of cyclic permutations was seen to be (p — l) n /n; only the fraction (1 — (3) n of 
these cannot be shortened by using the group identity, so we are down to (1 — (3 n ){p — l) n /n words. Roughly 
P n (k, a)(l — (3) n (p — 1)" /n of these, however, contain precisely k sequences of four letters each of which comes in two 
equivalent pairs due to the group identity; to avoid overcounting we have to divide out the multiplicity 2 fc . Finally 
summing over k we get the desired number of different n-letter words as 

Kn) = £ j^y^ii - - 0TiP - iT In 

= -[{l - a/2)(l - 0){p- 1)]" = I e " ln ^« , (B3) 
n n 

Peff = (l-a/2)(l-/3)(p-l) = 6.98. (B4) 

Interestin gly, the effective number of symbols is not drastically reduced from p = 8 — 1 = 7 by the group identity 
equation ( p.9| ). 
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